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O , Abstract 

o 

\Q • This is a gentle introduction to Colombeau nonlinear generalized func- 

tions, a generalization of the concept of distributions such that distributions 
^ | can freely be multiplied. It is intended to physicists and applied mathemati- 

^' cians who prefer a 'step-by-step approach' to a 'top-down indoctrination.' 

^ ■ No particular prerequisite knowledge is necessary — and in less than one 

hour you should know everything you need to know and were afraid to ask 
about Colombeau algebras and their applications in physics... 

A selected bibliography is appended, giving examples of applications 
pq . to partial differential and wave equations, electrodynamics, hydrodynamics, 

>> ' general relativity, and quantum field theory. 

o\ ■ 

(N 

in 

o 



The goal of this tutorial is to lead the reader to rediscover by himself the key 
ideas which led Colombeau to define the proper generalization of the concept of 



distributions such that multiplication is always possible and meaningful. 

00 : 

The emphasis is on concepts and methods, and the intent is to convince the 
reader that working with Colombeau nonlinear generalized functions (in short, 
(/-functions), which can be differentiated and multiplied freely, is not more com- 
S-h ' plicated than working with the familiar C°°-functions. 

a | 

Since everything is self-contained and kept simple there are only few references 
in the text. On the other hand, a selected bibliography with references to major 
publications on the subject is appended at the end. 

While Colombeau's seminal books 0, are still highly valuable, the most 
recent comprehensive textbook is |fTTfl. Short summaries of the main features of 
the Colombeau theory are included in most publications cited in the bibliography. 
Furthermore, an alternative primer on Colombeau algebras is given in lfT2ll . 



'The adjective 'nonlinear' emphasizes that Colombeau generalized functions form an algebra. 



1 Regular and irregular distributions 

The discovery of the Colombeau algebras is certainly one of the great events of the 
Twentieth Century history of mathematics. To understand how it came about, let 
us start from another great invention, that of the theory of distributions by Bochner, 
Sobolev, Mikusinski, and Schwartz. Indeed, whereas a regular distribution is a 
functional (f)(T) having the representation 

where (j)(x) is a locally integrable functional and T(x) a 'test' functionjj there is no 
such representation for the Dirac 'function' 5(x) which is defined by the functional 

5(T):=T(0), \/T(x)eV. (1.2) 

Thus, before the theory of such singular distributions was invented, the only thing 
that could be done was to write, symbolically, 



SW-JtoSW), VT W£ X>, ,,.3, 

and to refer to the definition (|1.2I) for the interpretation of (11.31) . 

2 The abstract and sequential views of distributions 

Schwartz showed that S(x) can be interpreted as an element of the space V of 
continuous linear functionals on V, and its derivatives defined as the derivatives 
of these functionals. That is, if 7 G V is any distribution, its derivatives in the 
'distributional sense' are such that, VT e V and D n = d n /dx n , 

/^D^wr^t-ir/^MD-rK*). (2,, 



2 In simple words, a function is locally integrable if it is integrable on every compact set. 

3 T> is the space X>($1) of C°° functions with compact support on an open subset ft C M. 
For simplicity of notation we write C, C m , and C°° for the continuous, m-times continuously 
differentiable, and respectively smooth functions with compact support on ft. We similarly write 
C p for the piece-wise continuous functions. Then C°° C C m C C° = C C C p . We also tacitly 
assume that all integrations are over K, and that all functions are extended to K by setting them 
equal to zero outside of Q. Finally, we set No = {0, N}. 



Alternatively, following Mikusinski, (11.31) can be written as the weak limit of a 
sequence of C°° functions S e , that is, 

5(T) := lim / dx S e (x)T{x) = T(0), VT(x) E V. (2.2) 

Indeed, if 5 e (x) is any family of functions 

6 e (x)=p e (x):=- e p{^), (2.3) 

where e e]0, 1 [ is a parameter, and p taken in the seflj 

A := \p(x) e S, and / dx p(x) = l|, (2.4) 

making the change of variable x = ey and taking the limit, it comes 

/ dx 5{x)T{x) := lim / dy p{y)T(ey) = T(0). (2.5) 

Returning to equation (12.31) one can observe that the sequence 5 e representing the 
^-distribution can actually by written as the convolution 



5 e (x) = I dy p e (x - y)5(y) = p e (x), (2.6) 

provided the symbol S inside the integral is interpreted according to its functional 
definition (fL2l). 



This method is general: It can be proved that convoluting a regular or singular 
distribution with any p e provides a representative sequence 

7 e (x) := p e (x) * 70) = dy p e (x - y)j(y), (2.7) 

of that distribution. As j e (x) E C°° this process of generating a smooth repre- 
sentative of j(x) E V is called a regularization, and the regularizing functions 
defined by (12.41) are termed regularizes or mollifiers. Thus, if 7 is any regular or 
singular distribution, (11.11) can be written 

7 (T) := lim / dx %(x) T(x), VT(x) E V, (2.8) 



4 S is the space of C°° functions with steep descent, i.e., such that f(x) £ S and its derivatives 
decrease more rapidly than any power of l/\x\ as x tends to infinity. In distribution theory one 
usually takes p G V c S because this enable to deal with distributions with non-compact support 
unrestrictedly. 



and 7 is then interpreted as the equivalence class of the weakly convergent se- 
quences of the smooth functions 7 e modulo weak zero-sequencies. 

In comparison to Schwartz's abstract theory, the advantages of Mikusinski's 
sequential view are that it provides explicit representations for the distributions, 
and that their derivatives are simply obtained by differentiating the representative 
sequencies. 



3 Schwartz's local structure theorem 

A particularly important contribution of Laurent Schwartz is the formulation of his 
local structure theorem stating that "any distribution is locally a partial derivative 
of a continuous function" [1, Theorems XXI and XXVI]: 

Theorem 1 Let T>'(Q) be the space of distributions on the compact set Q. Then 
every 7GD' is of the form 

1 (x) = J2v n 9n(x), (3.1) 

n 

where n E No, and the support of each g n E C(Q) is contained in an arbitrary 
compact neighborhood K C Q. 

For example, Dirac's function S(x) is generated by the second distributional deriva- 
tive of the absolute value \x\ E C°, i.e., S(x) = 1/2 D 2 |x|. 

Differentiation induces therefore the following remarkable cascade of rela- 
tionships: continuously differentiate functions — > continuous functions — > dis- 
tributions. This gives a unique position to Schwartz distributions because they 
constitute the smallest space in which all continuous functions can be differentiated 
any number of times. For this reason it is best to reserve the term 'distribution' to 
them, and to use the expression 'generalized function' for any of their generaliza- 
tions. On the other hand, classical generalizations of the concept of function such 
as piece-wise continuous functions, measures, Cauchy and Hadamard finite-parts 
of integrals, etc., are all distributions. 

For application in physics Schwartz's structure theorem is of great significance 
because it asserts that singular distributions do not come 'from nowhere,' but 
derive from a generating function g(x) E C°. For example, the classical electron 
charge distribution originates from the absolute value in the definition of the 



Coulomb potential, i.e., <\> = e/\r\, and due consideration to this fact leads to a 
distributionally consistent introduction of point charges and dipoles in classical 
electrodynamics ll25l |27l |28l l29l |30t 



4 Schwartz's multiplication impossibility theorem 

Distributions generalize ordinary functions, which can be regarded as trivial cases 
of distributions. They enjoy most of the properties of C°° functions (e.g., they can 
be differentiated any number of time) with the notable exception of multiplication. 

For example, if the product of distributions is defined in the most natural 
way, i.e., by multiplying representative sequencies, the square of the ^-function 
corresponds to (5 2 ) e (x) = (5 e ) 2 (x) = p 2 (x). Then, when evaluated on a test 
function T E V according to (12.81) . we get, 

f T(0) f 
dxS 2 (x)T(x) := lim / dx p 2 (x)T(x) = lim / dy p 2 (y) = oo, (4.1) 



which implies that 5 2 is not a distribution. Many mathematicians have of course 
tried to define a consistent product of distributions. But these efforts only con- 
firmed that there is no multiplication on all of V which still has values in V, 
unless some essential properties are given up. For instance, in any formulation 
such that the usual relations x ■ 1/x = 1 and x ■ 5(x) = hold, associativity leads 
to the contradiction 

(x-—)-5(x) = 5(x), whereas — ■ (x ■ 5(x)) = 0. (4.2) 

The goal therefore shifted towards finding an algebra Q of generalized functions 
containing the distributions and preserving most of the desirable properties of 
ordinary functions. But even that less ambitious goal turned out to be quite 
difficult. In particular, there are many options and it is not possible to know a 
priori which 'essential properties' should be preserved. For instance, possibly 
inspired by the cardinal position of continuous functions in his structure theorem 



(|3.1I) . Laurent Schwartz was particularly attached to the idea that these functions 
should have a similar position in Q. He therefore postulated a set of minimum 
requirements which can be phrased as follows^ 



Oj The differential algebra Q is associative and commutative. Its elements are 
written [u] when it is useful to emphasize that u e Q. 



5 See for example ifTTl p. 6]. 



1 1 1 The space of distributions V is linearly embedded into Q, and the function 
f(x) = 1 is the unit element for their product '0' in Q, i.e., V7 G V, there 

is an embedding V' — > (/, 7 1— >■ [7], and [1] [7] = [7]; 

_2_ There exists a derivation operator D : Q — ► C? that is linear and satisfies the 
Leibniz rule, i.e., Vw, v E Q, D(u u) = (Dm) © t> + u © (Du); 

[3] D restricted to [V] is the usual partial derivative consistent with the integra- 
tion by parts formula (12.11) ; 

[4] The product of two continuous functions embedded in £ coincides with the 
usual pointwise product '■' in C, i.e., V/, g E C, [f] [g] = / • g. 

Unfortunately, on the basis of simple counter-examples, it is easy to show that there 
is no associative and commutative differential algebra Q satisfying the require- 
ments \\\ - \A\. For example, they lead to the conclusion D 2 |x| = 0, whereas, as 
recalled above, D 2 |a;| = 25 (x) in distribution theory. This is the famous Schwartz 
impossibility theorem of 1954. 



5 Colombeau's breakthrough 

It was only in 1983 that Jean-Fran§ois Colombeau was able to show in a truly satis- 
factory manner that it is actually possible to construct associative and commutative 
algebras satisfying \Y\ - \3\ provided [4j is replaced byO 



4' I The product of two C°° functions embedded in Q coincides with the usual 
pointwise product '•' in C°°, i.e., V/, g G C°°, [/] Q[g} = f ■ g. 

Therefore, since C°° C C, it was by relaxing the requirement \4j that it became 
possible to move forward: As C°° functions have much more powerful properties 
than continuous functions in general, e.g., Taylor's theorem with remainder, the 
problem became manageable. 



6 The original discoveries of Colombeau were made in a different context, and arose from more 
abstract considerations. But their success can be traced to the emphasis given to C°° rather than to 
continuous functions in general, an emphasis which may have a deep physical significance. 



6 The embedding space 



Colombeau's axiom _4^_ combined with axioms \Y\ - [3] implies that Q contains 



C°° as a differential subalgebra: This opens the way to the possibility that Q could 
be similarly contained in a larger differential algebra £ such that its elements 
would be C°° in the variable x. Since the mollified sequencies j e (x) representing 
the distributions are precisely C°° in the variable x, this suggests to define £, the 
embedding space, as the set of mapsj3 

^:=((/l:Ax^l, 



{(/. 
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(r],x)^{f € ) e {r],x)j, (6.1) 



which are C°° functions in the variable x E f2 for any given Colombeau mollifier 
7] E A q , where A q C Ao remains to be specified, and which depend on the 
parameter e e]0, 1[ through the scaled mollifier 

r) e (x) := -Vy-J- (6-2) 

Obviously, £ is an associative and commutative differential algebra with unit 
(t], x) h-> 1 with respect to pointwise multiplication. It contains C°° as the subset 
of the maps (16.11) which do not depend on i], i.e., (f e ) £ (x) = /(x). 

The distributions / £ V are then embedded in £ as the convolutions^ 

(/ e ) e (x) := fc(-aO * /(x) = /" %(^)f(y) 

dz 7](z) f(x + ez), (6.3) 



where, in order to define Q C £, the Colombeau mollifiers i] E A q may need 
to have specific properties in addition to those implied by (12.41) . In particular, 
since C°° C V there are two distinct embeddings of C°° in £: Its embedding by 
(|6.3I) as a subset of V, and its direct inclusion by (16.11) according to the maps 



(fe)e(x) = f(x). To be consistent with axiom |_4j, the mollifiers rj E A q have 



thus to be such that [(/ £ ) £ ](x) = [f](x) = f(x) for all f E C c 



7 The notation (/ e ) e where < e < 1, which will be later abbreviated as f e , emphasizes that 
(f e )e is an element of £ rather than a usual representative sequence (12.7b . 

8 This definition due to Colombeau differs by a sign from the usual definition (12.71 ) of regular- 
ization. 



7 Embedding of C°° functions 

To find these additional properties we begin by studying the embeddings and 
products of C°° functions. We therefore calculate (16.31) for / G C°° and apply 
Taylor's theorem to obtain at once 

(f e ) e {x) = f{x)Jdzri{z) + ... (7.1) 

+ -J {n \x) I dz z n r 1 {z) + ... (7.2) 

71! J 

+ -A^yjdz z q+1 V (z) f(* +1 \x + &ez), (7.3) 

where p n \x) is the n-th derivative of f(x), and $ e]0, 1[. Then, since 77 6 <S and 
/ has a compact support, the integral in (17.31) is bounded so that the remainder is 
of order O x (e q+1 ) at any fixed point x. 

Moreover, if following Colombeau the mollifier 77 is chosen in the set 

A q '■— < T)(x) G Aq, and dz z n r](z) = 0, Vn = 1, ..., g k 

all the terms in (|7.2I) with n G [1, q] are zero and we are left with 

V/ G C°°, (/ e ) e (x) = /(x) + O^e^ 1 ). (7.5) 

Therefore, provided the set A q is not empty and q can take any value in N, it 
is possible to make the difference (f e ) e (x) — f(x) as small as we please even if 
e g]0, 1[ is kept finite. If we now consider a product of two C°° functions, it is 
easily seen that equation (17.51) immediately leads to 



(7.4) 



Vm,wgC°°, (u e ) e {x) ■ (v e ) e {x) =u(x)-v( y x) + O x (e q+1 ), (7.6) 

where the remainder O x (e q+1 ) is still as small as we please for any e g]0, 1[ if q is 
large enough. 

8 Colombeau mollifiers and Fourier transformation 

Colombeau proved that A q is not empty and provided a recursive algorithm for 
constructing the corresponding mollifiers for all q G N. He also showed [5, p. 7], 
[51 p.l 13], [33, p. 169] that due to the Fourier transformation identities 

dx r)(x) = t|(0), and / dx x n r](x) = (-i) n j-j;(0), (8-1) 



which are valid \/r] E S, the conditions (17.41) on the moments of rj(x) can be 
replaced by equivalent conditions on the derivatives of its Fourier transform rf(p) . 
Thus by taking for rj(x) any real functions such that rj(p) = 1 in a finite neighbor- 
hood of p = 0, one automatically satisfies the conditions (17.41) for any n e N, that 
is for g as large as we please. For this reason the set of mollifiers 

Ao := Ul(x) G S, suchthat rf{0) = l\. (8.2) 

is written Aoo- In this paper all Colombeau mollifiers will be taken in that set. 

For example, with rj e Aoo the Colombeau embeddings of any two polynomi- 
als, and the products of these embeddings, are identical to these polynomials and to 



their ordinary products. That is, axiom |4' | is identically satisfied for polynomials. 



But for the other C°° functions there will still be a remainder to be taken care of, 
even if it is infinitesimal. 



9 Embedding of continuous functions 

Let us now consider the embedding of continuous functions assuming that the only 
things that are known is that they are continuous and compactly supported. Then, 
a priori, there is little more that can be done than writing 

V/ e C, (/ e ) e (z) = J dz n(z) f(x + ez), (9.1) 

because neither Taylor's formula nor the mean-value theorem can be applied to 
transform the right-hand side into a more useful expression. In fact, the only fully 
general expression comparable to (|7.5I) is 



V/eC, (/«)«(*) =/(z)+o* l6 (l), (9.2) 

which simply means that (f e ) e (x) converges uniformly to f(x) as e — > because 
/ has compact support. Any more precise statement requires that the continuous 
function is further specified. 

For example, if / is m-times continuously differentiable we can write 

V/eC"\ {f e Ux) = f{x) + O x {e m ), (9.3) 

where, in contrast to (17.51) . m > 1 is a fixed integer. 



10 Embedding of distributions 

We now turn to distributions. By Schwartz's local structure theorem (|3.1I) they 
can be written 7 (x) = D n g(x) where g E C(K) if we restrict ourselves to a 
single generating function with support in a compact set K. Then for x E K their 
embeddings (16.31) are, using the integration by parts formula (|2.1I) . 



( 7e ) e (x) = J dz rj(z) B n x g(x + ez), (10.1) 

= ±Jdz V (z)D n z g(x + ez), (10.2) 

= (— Y J dz (D n r?) (z) g[x + ez). (10.3) 

Since r\ E S, and g E C is compactly supported, the last integral is bounded and 
we get 

V 7 Gl?', ( 7e ) e (x) = !B (l/e n ), as e - 0. (10.4) 

This bound is compatible with the bounds (17.51 I9.2L and 19.31) because V contains 
all continuous functions. To illustrate its significance for non-trivial distributions 
we need to consider generating functions g E C°. 

For example, we know that S(x) = D 2 g{x) with g(x) = \x\/2. On the other 
hand, the Colombeau embedding (16.31) of S(x) is 

(Wx) = ^(~)=0„(l/e), (10.5) 

which has an e-dependent bound consistent with the bound (110.41) . i.e.. 1/e < 1/e 2 , 
although their exponents disagree by one unit. This is because (110.41) is fully 
general and thus does not take the particular properties of g{x) into account. In 
the present case it is easy to calculate (|a;| e ) £ (a;) with (16.31) and to verify that 
2(S 6 ) e (x) = D 2 (\x\ t ) t (x) = O x .(l/e) rather than O x (l/e 2 ). For the same reason 
the embedding of the Heaviside function H(x) = D|x|/2 is 

(H e ) t (x) = [ X( dz ri(-z) = 0,(1). (10.6) 



which, rather than O z (l/e), has the e-dependence 0^(1) characteristic of a piece- 
wise continuous function because of its jump at x = 0. 

To give another example, the singular distributions generated by the derivatives 
of the C° function equal to for x < and to x r for x > 0, with r e]0, 1 [, have an 
e-dependent bound O x (e r ~ n ) where r — n < in BL 

10 



Local structure of embeddings and of their differences in 6 


f 


(fe)e(x) 




c°° 


f(x) 


Directly included smooth function 


c°° 


f(x) + O x {&) 


Smooth function, q > Vp G N 


c 


f(x)+o x>e (l) 


Continuous function 


c 


0,(1) 


Piece-wise continuous function 


V 


O x (e~ N ) 


Singular distribution, N E N 


M 


O x (e") 


Negligible function, q > Vp G N 



Table 1 : The differential algebra £ contains the smooth functions as direct 
embeddings f G C°° C £, and also as mollified embeddings (/ 6 ) e G (Cf D ) e C £. 
The embeddings (/ 6 ) e o/f/ie continuous functions and of the distributions are sorted 
in terms of the behavior of the bound on their e-dependent part as e — > 0. /(x) « 
the point-value of the continuous functions at e = 0. 77ze negligible functions are 
infinite simally small elements such as the differences between the direct inclusions 
and the Colombeau-mollified embeddings of the C°° functions. 



In summary, (110.41) provides a conservative bound for the e-dependence of all 
Schwartz distributions. In the case of singular distributions, the exponent in the 
bound (110.41) can be any integer n G N. 



11 Linear operations and negligible functions 

In Tabled] the usual functions and the distributions are classified according to the 
structure of their embeddings in 8. Referring to this table it is easy to predict the 
structure of the result of binary algebraic operations in £, and thus to get clues on 
how to define the algebra Q. 

For instance, in the last line of Table Q] the set denoted by J\f consists of 
functions which are not the direct result of embeddings: It is the algebra of the 
so-called negligible functions^ which arise in particular from subtracting the two 
different inclusions of the C°° functions, i.e., 



V/GC°°, VgGN, (feUx)-f(x) = O x ((*) GAT. 



(11.1) 



A proper definition of negligible functions will be given shortly. 
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Then for all linear operations (i.e., addition/subtraction, multiplication by a scalar, 
and differentiation) it is clear that the results will always be in one of the sets 
listed in the first column of Table [H which is therefore a suitable classification 
of the embeddings of the usual functions and distributions with regards to linear 
operations in £. 

Of course, the negligible functions of the type (111.11) are precisely the differ- 



ences that are to be taken care of in order to satisfy axiom |4' | . In particular, 
they will remain 'negligible' as long as they are not multiplied by 'very large' 
functions. This is why we have now to look at the nonlinear operations in £. 



12 Nonlinear operations and moderate functions 

We know that the product of two distributions (or of a continuous function and a 
distribution) will, in general, not be a distribution. For example, the n-th power 
of the Dirac (5-function can be defined by the n-th power of the embedding (110.51) . 
i.e., 

(Oe(s) = ^ B (-~)=0.(e-*). (12.1) 

But, despite that (6™) e (x) has a O a ,(e _n ) dependence similar to that of a 'distribu- 
tion,' it is not a distribution in the sense of Schwartz and Mikusinski — rather, it 
is an element of a larger set of 'generalized functions' containing the distributions 
as a subspace. 

This led Colombeau to define the set £ M , which he called moderate functions)^ 

V(g e ) e E £ M : 3iV G N„, such that (g e ) e (x) = O x (e' N ). (12.2) 

It is evident that M C (C°°) e C (C) e C (D') e C £ M , and a matter of elementary 
calculations to verify that £ M and J\f are algebras for the usual pointwise operations 
in £, and that J\f is an ideal of £ M . Indeed, the product of two moderate functions is 
still moderate — they are multipliable — and as q in (|ll.ll) is as large as we please, 
and N in (112.21) a fixed integer, the product of a negligible function by a moderate 
one will always be a negligible function. Moreover, £ M is a differential algebra 
satisfying axioms \Y\ - \3J, and it is not difficult to show that £m is the largest 
differential subalgebra (i.e., stable under partial differentiation) of £ in which N 
is a differential ideal. 



'A proper definition of moderate functions will be given shortly. 
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Multiplication in £ 




c°° 


Af 


£m 


£m\£ 


c°° 


c°° 


Af 


£m 


£m\£ 


Af 


M 


Af 


M 


£ 


£u 


£m 


M 


£m 


£ 


£m\£ 


£m\£ 


£ 


£ 


£m\£ 



Table 2: The elements of £ remain in their respectives subalgebras Af, £m, or 
£ M \£ when multiplied by directly included C°° functions. The subalgebra Af C £m 
is an ideal of£f4- The products of negligible and moderate elements with elements 
in the complement of £m in £ are in general undefined. 



Furthermore, one can also consider infinite sums of products of moderate 
functions and take their limits in £ . It is then easily verified that, for example, 
\Z(SfY e (x) and sm(5™) e (x) are elements of £ M . On the other hand 



exp(|(J| e ) 6 (0) = O !B (exp(l/e))^£ ] 



M- 



(12.3) 



so that exp(S e ) e (x) is a non-moderate function, and thus an element of the com- 
plement of £ M in £. Conversely, exp(—\5\ e ) € (x) is a negligible function, so that 
Af contains elements of exponentially fast decrease. 

Consequently, when operating in full generality in £ , that is when including 
multiplication and limiting processes, one is led to consider its elements as in 
Table El i.e., as members of the differential algebras C°°, A/", and £ M , rather than 
as members of the embeddings of the classical spaces C°°, C, and V as in Table[TJ 



13 Discovery of the Colombeau algebra 

The fact that Af is an ideal of £ M i s th e key to defining an algebra containing the 
distributions and satisfying axiom _4^_. Indeed, if we conventionally write Af for 
any negligible function, then_| 



Vg e , K e E M , (g e + Af) ■ (h e + Af) = g e ■ h e + Af. 



(13.1) 



'From now on we abbreviate (/ e ) e as f e . 
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Similarly, using the same convention, Eq. (|7.6I) giving the product of the Colombeau 
embeddings of two C°° functions can be written 

Vu,veC°° } [u + M)-{v + M)=u-v+M, (13.2) 



whereas axiom |4' | demands 

Wu,veC°°, [u + N]o[v+Jtf\ = [u-v + JV]=u-v. (13.3) 

Thus, it suffice to define the elements of Q as the elements of £m modulo TV, e.g., 
to identify [g e + TV] and [g e ], so that [u + J\f] = [u] because (C°°) e C £m, and 
axiom 



will be satisfied. 



This immediately leads to the definition of the Colombeau algebra as the 
quotient 

Q ■= ^ (13-4) 

That is, an element g E Q is an equivalence class [g] = [g e + J\f] of an element 
g e G £m, which is called a representative of the generalized function g. The product 
g h is defined as the class of g e ■ h e where g e and h e are (arbitrary) representatives 
of g and h; similarly Dg is the class of Dg e if D is any partial differentiation 
operator. Therefore, when working in Q, all algebraic and differential operations 
(as well as composition of functions, etc.) are performed component-wise at the 
level of the representatives g e . 

Q is an associative and commutative differential algebra because both £ M and 
J\f are such. The two main ingredients which led to its definition are the primacy 
given to C°° functions, and the use of the Colombeau mollifiers for the embeddings. 



14 Special and general Colombeau algebras 

Depending on the precise definitions of the moderate and negligible functions, 
as well as of any further specification constraining the Colombeau mollifiers, 
there can be many variants of Q, even if the domain and range of the generalized 
functions are simply a subset of R. There are however two general types of 
Colombeau algebras: The 'special' (or 'simple') algebras, and the 'general' (also 
called 'full' or 'elementary') algebras. 

For example, let us define a special Colombeau algebra of generalized func- 
tions on fi 6 R n with value in C. Then, using the standard multi-index notation 

0|«I 

D a = (14 1) 

{dx x )^---(dx n )<*^ K ' ; 
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where aGNJ and |a| = ai + a 2 + • ■ ■ a n , the distributions will be the partial 
derivatives D a f(x) of the continuous function f{x) G C{Vt). A possible definition 
of G s (Vt), which can easily be adapted to more complicated manifolds, is as 
follows: 

Definition 1 (Embedding space) Let Vt be an open set in R n , let e g]0, 1[ be a 
parameter, and let 77 G Aoo be an arbitrary but fixed Colombeau mollifier. The 
'embedding space ' is the differential algebra 

£ 8 (ty:={f e : ]0,l[xfi-C, 

(e,£)h-/ e (a?)}, (14.2) 

where the sequencies f^are C°° functions in the variable x G Q. The compactly 
supported distributions^ are embedded in £ s by convolution with the scaled mol- 
lifier r] £ , i.e., 

f e (x) := J^fv ( IL jp) f(y) = J dz n V (z) f(x + ez). (14.3) 

Definition 2 (Moderate functions) The differential subalgebra £^ C £ s of 'mod- 
erate functions' is 



£^(0) := <f e - VK compact in Q, Va G N£, 
3N G N swc/z ?/za?, 
sup \D a f e (x )| = 0(4f) a* e -► O). (14.4) 

Definition 3 (Negligible functions) The differential ideal M s C £^ of 'negligible 
functions ' is 

J\f*(Q) := |/ £ : Vif compact in Q, Va G N™, 
VgGN, 
sup |D a / e (f)| = 0(e 9 ) ase^o). (14.5) 

Definition 4 (Special algebra) The special Colombeau algebra is the quotient 

FPU := $g. (14.6) 



12 The embedding of all of V is achieved by a more complicated formula based on sheaf theoretic 
arguments, see [11, Proposition 1.2.13] . 
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The main differences with the 'naive' definitions (111.11) and (112.21) are: (i) The 
resort to the compact subset K C Vt, which is necessary because of the local 
character of the concept of distribution, as is clearly stipulated by Schwartz's 
structure theorem; (ii) the need to consider the supremum over all K C Q in 
order to take into account all possible discontinuities when x ranges in fi; and (iii) 
the need to consider all possible derivatives of f € in order that the moderate and 
negligible functions have the required properties for all their derivatives. 

A general Colombeau algebra Q g is an enlargement of Q s , obtained by consid- 
ering all r] E Aqo and by replacing (in both J\f s and E^) the functions x — * f e (x) 
by the set of functions x — > f e (x, rf) depending on r]. Since all possible 77 are con- 
sidered the arbitrariness characteristic of £ s (f2) disappears, and the embeddings 
of the distributions and functions with finite differentiability become 'canonical' 
since they do not depend any more on a fixed mollifier. However, while this is 
conceptually interesting from the mathematical point of view, it is not a necessity 
since the particular mollifier (or set of mollifiers) defining a special Colombeau al- 
gebra Q s C Q 9 may have a physical interpretation. For this reason the dependence 
of the embeddings on the mollifiers is not a defect, but rather a positive feature in 
many applications of the special Colombeau algebras. 



15 Interpretation of distributions within Q 

To construct the Colombeau algebra we have been led to embed the distributions 
as the representative sequences (7 e ) 6 E £ defined by (16.31) where r\ e E Aoo is a 
Colombeau mollifier, that is not as the usual representative sequencies defined 
by (|2.7I) where p e E A . However, since A^ C Aq, we can still recover any 
distribution 7 from 7 e = (j e ) e by means of (12.81) . i.e., as the equivalence class 

7 (T) := lim / dx j e (x) T(x), VT(x) E V, (15.1) 

where 7 e can be any representative of the class [7] = [7^ + TV] because negligible 
elements are zero in the limit e — >• 0. 

Of course, as we work in Q and its elements get algebraically combined with 
other elements, there can be generalized functions [g e ] different from the class [%] 
of an embedded distribution which nevertheless correspond to the same distribution 
7. This leads to the concept of association: We say that two generalized functions 
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g and h are associated, and we write g x hj^j iff 

lim J dx (g e (x) - h e (x)\T(x) = 0, VT{x) e V. (15.2) 

Thus, if 7 is a distribution and g some generalized function, the relation 9x7 
implies that g admits 7 as 'associated distribution,' and 7 is called the 'distributional 
shadow' (or 'distributional projection') of g because the mapping 7 e h-> 7 defined 
by (115.11) is then a projection C? — > £>' for all g e associated to j e . 

Objects (functions, numbers, etc.) which are equivalent to zero in Q, i.e., equal 
to 0(e 9 ), Vg G N, are called zero. On the other hand, objects associated to zero in 
Q, that is which tend to zero as e — > 0, are called infinitesimals. Definition (115.21) 
therefore means that two different generalized functions associated to the same 
distribution differ by an infinitesimal. 

16 Multiplication of distributions in Q 

The continuous functions and their derivatives, i.e., the distributions, are not 
subalgebras of Q: Only the smooth functions have that property. Thus we do not 
normally expect that their products in Q will be associated to some continuous 
functions or distributions: In general these products will be genuine generalized 
functions, i.e., new mathematical objects — which constitute one of the main 
attractions of Q. 

For example, the n-th power of Dirac's (^-function in Q, Eq. (112.11) . has no 
associated distribution. But 5 n is a moderate function and thus makes perfectly 
sense in Q. Moreover, its point- value at zero, rj n (0)/e n can be considered as a 
'generalized number.' 

On the other hand, we have elements like the n-th power of Heaviside's 
function, Eq. (110.61) . which has an associated distribution but is such that [H n ] (x) 7^ 
[H](x) in Q, whereas H n (x) = H(x) as a distribution in V . Similarly, we have 
[x] [5](x) 7^ in Q, whereas xS(x) = in V. In both cases everything is 
consistent: Using (115.21) . one easily verifies that indeed [H n ](x) x [H](x) and 
[x}Q[5](x) xO. 

These differences between products in Q and in V stem from the fact that 
distributions embedded and multiplied in Q carry along with them infinitesimal 



3 In the literature the symbol w is generally used for association. We prefer to use x because 
association is not some kind of an 'approximate' relationship, but rather the precise statement that 
a generalized function corresponds to a distribution. 
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information on their 'microscopic structure.' That information is necessary in 
order that the products and their derivatives are well defined in Q, and is lost when 
the factors are identified with their distributional projection in V . 

Let us illustrate this essential point with a concrete example. In physics the 
Heaviside function H(:r) represents a function whose values jump from to 1 
in a tiny interval of width e around x = 0. Thus it is obvious that J(H 2 (x) — 
H(x))T(x) dx tends to when e — > + if T is a bounded function, i.e. ,H 2 x H. But 
since H' is unbounded one has f (H 2 (x) - H(x)) .H'(x) dx = 1/3 - 1/2 = -1/6, 
as obvious from elementary calculations. This shows that one is not allowed to 
state H 2 = H in a context where the function H 2 — H could be multiplied by a 
function taking infinite values such as the Dirac function S = H'. 

Therefore, the distinction between ^-functions that are 'infmitesimaily nonzero' 
such as H 2 — H from the genuine zero function insures that multiplication is co- 
herent in Q, because 'infmitesimaily nonzero quantities,' when multiplied by 
'infinitely large quantities,' can give significant nonzero results. At the same time, 
this distinction insures that all calculations are consistent with those in V . In 
particular, if at any point it is desirable to look at the intermediate results of a 
calculation from the point of view of distribution theory, one can always use the 
concept of association to retrieve their distributional content. In fact, this is facili- 
tated by a few simple formulas which easily derive from the definition (115.21) . For 
instance, 

V/i,V/ 2 eC =► [/i]0[/ 2 ]x[/rg (16.1) 

and, 

V/gC°°,V 7 gX?' => [/]0[t]x[/- 7 ], (16.2) 

but, in general, 

V71, V72 e V => [71] [ 72 ] t [71 • 72], (16.3) 

whereas, 

V<?i,V# 2 6£, 91-92 => D ft Sl xD a j, (16.4) 

For example, applying the last equation to [H 2 ](:r) x [H](x) one proves the often 
used distributional identity 2 [5] (x) [H] (x) x [5] (x) . 

In summary, one calculates in Q as in C°° by operating on the representatives 
g e £ £ with the usual operations {+, — , x, d/dx}. The distributional aspects, if 
required, can be retrieved at all stages by means of association. 
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17 Working with distributions versus working with 
^-functions 

An interesting feature of Colombeau algebras is that they enable, in many cases, 
to set aside the concept of distributions and to replace it by the more general and 
flexible one of ^-functions. 

Indeed, a distribution cannot be the end result of a calculation in any physical 
theory: It is a functional which has to be integrated over its argument to yield 
a quantity comparable to experiment. Similarly, a measurement is always some 
kind of an average over a continuous distribution of matter supported by bodies of 
finite extension. Thus, if one takes the sequential view, one is often led in physics 
to consider integrals of the type (12.81) . i.e., 



g(S) = lim f dx g e (x) S(x), (17.1) 

where g(x) may be any regular or singular distribution corresponding to a basic 
physical quantity (e.g., an energy density), and S(x) E V a smooth function (e.g., 
S(x)dx could be a volume element). 

There are then two options: 

• In conventional 'distribution theory' the distributional aspect is emphasized 
throughout the calculation and all intermediate results are interpreted as 
distributions. In particular, when working according to the sequential view, 
limits similar to that in (117.11) are taken at all stages so that information 
that could be relevant to nonlinear operations is discarded. (In the language 
of generalized functions, one systematically works with the distributional 
shadows rather than with the generalized functions themselves.) The method 
is therefore restricted to linear theories, and if the limit e — > is undefined 
the end result will in general be meaningless even if e is kept finite. (Because 
infinitesimal information that could have been significant before passing to 
the limit may have been discarded.) 

Example: If the electrostatic Coulomb potential is defined as a distribution, 
the ^-embedded Coulomb field has the form ||24ll27ll30H 

EJr ) = e limf— H(r — a)) — elimf— -5(r — a)) (17.2) 



r 



elim— H(r-o). (17.3) 
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In distribution theory, the second part of (117.21) . which contains the 5- 
function, is discarded 11241 p. 144]. Expression (117.31) cannot therefore be 
used in non-linear calculations such as the self-energy of an electron. 

• In '(/-function theory' all non-smooth functions / are represented by their 
Colombeau mollified sequence f £ , and there is a unique parameter e which 
is kept finite until the end of the calculation. (There is also possibly a unique 
common mollifier 77 if one works in a special Colombeau algebra.) All 
operations are then performed on these representatives, and at any stage one 
can verify the validity of the calculations by checking that the intermediate 
results are moderate functions. It is possible to consistently manipulate sin- 
gular distributions in nonlinear calculations, and the end results are obtained 
by taking the limit e — > as in (117.11) . If the theory is linear, these results 
are identical to those of the conventional 'distribution theory' option. In 
linear or nonlinear theories which lead to divergent quantities as e — ► the 
parameter e can be left finite, and the end result can be interpreted as a 
'generalized number.' This generalized number may then be renormalized 
to some finite quantity, which implies that any dependence on e and on the 
arbitrary mollifier i] is removed at this final stage. 

Example: Calculating the self-energy of an electron involves integrating 
the square of its electric field. In distribution theory, where this square is 
undefined, the result using (117.31) is the well-known expression 

U seif (V') = ?-\im-, (17.4) 

Z a^O a 

which diverges as a — > 0, and which corresponds to the square of the first 
term in (117.21) . i.e., to the energy of the field surrounding the electron. On 
the other hand, when calculated in Q using (117.21) the self-energy is ll27l[30l 

2 1 /•+£» 

U seli (g) = -\im- dx V 2 (-x). (17.5) 

2 e->o e J_ OQ 

which is independent of the 'cut-off' a, and which corresponds to the square 
of the second term in (117.21) . i.e., to the square of a ^-function. Thus, when 
calculated in Q, the self-energy is entirely located at the position of the 
electron, i.e, precisely where its 'inertia' as a point-mass is supposed to 
reside. It remains therefore to renormalize U se \f(Q) to the measured mass of 
the electron, and everything makes mathematically and physically sense. 

That discussion permits to conclude this paper by an analogy: The relations of 
the usual functions and distributions to the (/-functions are somewhat analogous 
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to those of the real to the complex numbers. If one looks at e as an analog 
of i, then taking the real part of a complex number corresponds to restricting a 
generalized function to its associated function or distribution by taking the limit 
e — ► 0. Working in M. or V is therefore less general and flexible than working in 
Cor£. 



Acknowledgments 

The author is very much indebted to Drs. G. Falquet and J. -P. Hurni for stimulating 
discussions and comments and for a critical reading of the manuscript. 



21 



Appendix: Special algebra of tempered £/ -functions 

In many mathematical and physical applications the restriction of the test functions 
to the space V of C°° functions with compact support is too restrictive 00. This 
is the case of the Fourier transform which, in its simplest form, has as a kernel 
cos(px) which is not integrable over the whole space fi = R. Thus, just like in 
Schwartz distribution theory, an extension of the Colombeau theory to 'tempered 
(/-functions,' see, e.g., IfTTl p. 15 and p. 65], is essential when dealing with func- 
tions that are integrated over the whole space fi = R n . Moreover, the algebra Q l 
of tempered (/-functions has a property that is important from a practical point of 
view: In Q l , componentwise composition is a well defined operation generalizing 
composition of C°° functions. 

To use this extension it is necessary to be careful about the definitions of the 
pertinent function spaces. We therefore recall IfTTl p. 15]: 

Definition 5 (Algebras S, O c , and O m ) Let Q l c R n be a n- dimensional box, 
i.e., a subset of the form iO\ x • - • x u n where Ui is a finite or infinite open interval 
in R. Then, 

S ■= If G C°°(fi 4 ) : Vm G N,Va G N£, 

sup (l + |f|) +m |D a /(f)| <oo), (18.1) 

Sen* ' 



O c := {/ e C°°(^*) : 3m G N such that, Va e N™, 

sup (l + |f |)- m |D a /(f)| <oo), (18.2) 

sen* J 

O m := {/ G C°°(ft*) : Va G N™, 3m G N, 

sup (l + |f |)- m |D Q /(f)| <oo). (18.3) 

sen* J 

Oq and O m correspond to two closely related definitions of functions with poly- 
nomial growth as \x\ — > oo. But, while O m corresponds to the usual definition, it 
is the algebra O c which in the (/-context provides the proper 'tempered' extension 
of the notion of C°° functions with compact support. 

Definition 6 (Embedding space of temperate distributions) Let Q l C R n be a 

n-dimensional box, let e g]0, 1[ be a parameter, and let r\ G Aoo C S be an 
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arbitrary but fixed Colombeau mollifier. The 'embedding space ' is the differential 
algebra 

£<(&):={&: ]0,l[xtf-C, 

(e,£)i-/ e (£)}, (18.4) 

where the sequencies f e are C°° functions in the variable x G Vt l . The tempered 
distributions g G S' are embedded in £ t by convolution with the scaled mollifier 
r] e , i.e., 

g e (x) := j^f v(^) 9(y) = j dz n r)(z) g{x + ez). (18.5) 

The functions h G Oq are directly embedded, i.e., 

V/i G O c , (h e ) e (x) = h(x), (18.6) 

so that O c is a subalgebra ofS 1 ^). 

Definition 7 (Temperate moderate functions) The differential subalgebra 8^ C 
£ f of 'moderate functions' is 

£^(ft*) := |/ e : VK compact in to*, Va G N£, 
3iV G N swc/i that, 
sup (l + \x\)- N \D a f e (x)\=0(-^) ase^o). (18.7) 

Definition 8 (Temperate negligible functions) The differential ideal J\f l C S^ 
of 'negligible functions' is 

A/"'(^) := |/ e : VfsT compact in to*, Vet G N \ 

3m G N swc/j ^a?, Vg G N, 

sup (l + |f |)~ m |D a / e (f)| = 0(e 9 ) ase^o). (18.8) 
xeK > 

Definition 9 (Special algebra of tempered ^-functions) The special Colombeau 
algebra is the quotient 

OW - ||}. (18.9) 
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